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ABSTRACT 

Stellar dynamical modeling is a powerful method to determine the mass of black holes in quiescent galaxies. 
However, in previous work the presence of a dark matter halo has been ignored in the modeling. Gebhardt & 
Thomas (2009) showed that accounting for a dark matter halo increased the black hole mass of the massive 
galaxy M87 by a factor of two. We used a sample of 12 galaxies to investigate the effect of accounting for a dark 
matter halo in the dynamical modeling in more detail, and also updated the masses using improved modeling. 
The sample of galaxies possesses Hubble Space Telescope and ground-based observations of stellar kinematics. 
Their black hole masses have been presented before, but without including a dark matter halo in the models. 
Without a dark halo, we find a mean increase in the estimated mass of 1.5 for the whole sample compared 
to previous results. We attribute this change to using a more complete orbit library. When we include a dark 
matter halo, along with the updated models, we find an additional increase in black hole mass by a factor of 
1.2 in the mean, much less than for M87. We attribute the smaller discrepancy in black hole mass to using data 
that better resolve the black hole's sphere of influence. We redetermined the M. - <x, and M, - Ly relationships 
using our updated black hole masses and found a slight increase in both normalization and intrinsic scatter. 
Subject headings: black hole physics - galaxies: general - galaxies: kinematics and dynamics - galaxies: nuclei 



1. INTRODUCTION 

It is now well established that almost every massive 
galaxy harbors a supermassive black hole in its center. 
Furthermore, close relations between the mass of this su- 
permassive black hole and the properties of the galaxy's 
spher oid component have been found, namely with the 
mass dMagorri an et al.l 1 1 998t lHaring & Rixll2004|). luminos- 
ity ([Kormendv & Richstong 119951; iKormendv & Gebhardll 
200 ll iMarconi & Hunt! I2003t iGultekin et alJ l2009bl> 
and w ith the velocity disp e rsion ^Gebhardt et all 
2000al: iFerrarese & Merrittl 12000: iTremaine et alJ 120021; 



Giiltekin et al.ll2009bl) . These relations imply a link between 
the growth of black holes and galaxy evolution, usually 
attributed to active galactic nucleus ( AGN) feedback (e.g 



Silk & Reesl 119981: bi Matteo et alJ 120051: ISpringel et ail 



20051: ICiotti & Ostrikerf 12007^ but at least to some degree 



they are a natural r esult within a merger-driven galaxy 
evolution framework (|Pend 120071: iHirschmann et al] 12010c 
iJahnke & Macciol 120101) . In general, the black hole-bulge 
relations and especially their evolution with cosmic time are 
able to provide deep insight into galaxy formation and black 
hole growth. Therefore, it is essential to properly establish 
the local relationships as precisely as possible. 

The black hole-bulge relations are based on a sample of 
~ 50 quiescent black holes, whose m asses have been deter- 
mined based on maser emission (e.g., iGreenhill et all 120031: 
Herrnstein et al. 2005; Kuo et al. 20 10[). gas kinematics (e.g., 
Ferrarese et al. mT9M Kconi et al.l 1200 lh iDalla Bonta et al. 



20091). and stel l ar kinematics (e.g., van der Marel et al.ll 19981; 
Gebhardt et all 120001* IShapiro et al.l 120061; iGebhardt et al] 
20071: IGultekin et al] l2009ah . In particular, stellar dynami- 
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cal modeling using orbit superposition is a powerful method 
to estimate black hole masses in quiescent galaxies, also re- 
covering the orbital structure within the galaxy. Usually, ax- 
isymmetry is assumed in these models. However, there are 
still uncertainties and possibly systematic biases within these 
methods. Uncertainties may arise from the deprojection of the 
observables onto three dimensions as the true inclination of 
the galaxy often is not well known, the presence of dust, some 
triaxiali ty that cannot be modeled properly with axisymmetric 
models (Ivan den Bosch & de Zeeuwl2010l) . or the presence of 
an AGN at the center. So far, in most models the contribution 
of the galaxy's dark matter (DM) halo has been neglected. 
IGebhardt & Thomasl ( 120091) recently showed that the black 
hole mass can be underestimated in this case. For the mas- 
sive galaxy M87 they found an increase of more than a factor 
of two in the measured black hole mass, just by including a 
DM halo in the modeling. The reason is that without a DM 
halo the mass-to-light ratio is overestimated in order to ac- 
count for the mass in the outer parts of the galaxy. Under the 
usually applied assumption of a constant mass-to-light ratio 
for the whole galaxy, this will propagate inward and lead to 
an underestimation of the black hole mass at the center due to 
overestimation of the stellar contribution. 

A s imilar result has been obtained by iMcConnell et al] 
( 120101) . They measured the black hole mass in the brightest 
cluster galaxy NGC 6086. They report a factor of six dif- 
ference between the black hole mass obtained from models 
without a DM halo and their most massive DM halo mod- 
els. However, the black hole's sphere of influence is barely 
re solved in their work. 

IShen & Gebhardtl (120101) found for NGC 4649, also a mas- 
sive galaxy, no change in the black hole mass by including a 
DM halo. In this case the sphere of influence is well resolved 
by the data. A larger sample, especially spanning a larger 
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ra nge in mass, is c learly required. 

IGebhardt et alj (12003b (hereafter G03) studied a sample of 
12 galaxies with kinematics derived from Hubble Space Tele- 
scope (HST) and ground-based observations, using axisym- 
metric orbit superposition models. They do not include a DM 
halo in their modeling. Since then the orb it superposition code 
used by our group has been improved dThomas et al.ll2004l 
l200HlSiopisetal.ll2009l) . by including a more complete orbit 
sampling. 

The aim of this paper is to reanalyze the data set presented 
by G03, using the most up-to-date dynamical modeling code 
and investigating the effect of accounting for the dark mat- 
ter contribution on the derived black hole masses. One of the 
galaxies in the G03 sample, NG C 4649, was recently ana- 
lyzed bv lShen & Gebhardt! (120101) , including a DM halo. We 
have reanalyzed this galaxy for consistency with the remain- 
ing sample, but find consistent results with this previous in- 
vestigation. 

2. DATA 

The data set used in this work is identical to those used in 
the work of G03. Thus, we will only give a brief summary 
and refer to G03 for more detail. The data consist of three 
sets of observations for each galaxy: imaging, HST stellar 
kinematics, and ground-based stellar kinematics. 

High-resolution imaging is required to obtain the stellar sur- 
face brightness profile for each galaxy. This imagi ng has been 
obtain ed in the V band with the HST WFPC2 (lLauer et alj 
[20051) . except for NGC 4697, w hich was observed with the 
HST WFPC1 (lLauer et al.ll 19951) . Surface brightness profiles 
were measured from the pont-spread function deconvolved 
images and were augmented with ground-based imaging at 
the outer parts, not covered by HST. For the deprojection of 
the surface brightness profile to a luminosity density profile, 
we assume axisymmetry, an inclination angle of 90°, which 
we refer to as edge-o n, and used the technique outlined in 
IGebhardt et all (fl996l) . 

The HST observa tions and kinematics are presented by 
iPinknev et alj ([2003) and G03. They consist of Space Tele- 
scope Imaging Spectrograph (STIS) long-slit spectra along 
the major axes, except for NGC 3377 and NGC 5845, 
which have Faint Object Spectrograph (FOS) aperture spec- 
tra. The spectra cover the Ca II triplet around 8500 A. For 
the dynamical modeling, the line-of-sight velocity distribu- 
tions (LOSVDs), extracted from the spectra, are used. The 
LOSVDs are given in a non-parametric form, binned into 15 
equidistant bins, compared to 13 bins in G03. 

The grou n d-base d kinematics are presented by 
IPinknev et al.1 ( 120031) and G03 as well. They consist of 
long-slit spectra along several position angles, mainly ob- 
tained at the MDM observatory. They also cover the Ca II 
triplet, or alternatively the Mgb absorption at 5175 A. The 
individual LOSVDs are binned into 15 points as well. 

3. DYNAMICAL MODELS 

The dynamical modeling uses the orb it supe rposition 
method, first proposed by ISchwarzschildl d 1979b . This 
general m ethod h as been widely used by v arious groups 
(iRix et al.1 119971: Ivan der Marel et al.1 11998b ICretton et all 
ll999tlValluri et al.ll2004 ). Our techn ique is described i n detai l 
in G03, iThomas et alj O004l[2005h and lSiopis et ail (120091) . 
We will give a brief summary here and especially point out the 
differences compared to the work of G03. The basic approach 
consists of the following steps: (1) deprojection of the surface 



brightness profile to a three-dimensional luminosity distribu- 
tion, (2) computation of the specified gravitational potential, 
(3) generation of a representative orbit library in this poten- 
tial, (4) fitting the orbit library to the observed light distribu- 
tion and kinematics, and (5) modifying the input potential to 
find the best match to the data, based on &x~ analysis. 

As described in Sec tion [2] we depro j ect th e surface bright- 
ness profile following IGebhardt et al.1 (| 19961) and assume an 
edge-on configuration, as used by G03. The only exception 
is NCG 4473, where we assume an inclination of 72°, as has 
been done in G03. To determine the potential, we assume a 
constant mass-to-light ratio throughout the galaxy, a specific 
black hole mass and optionally also include a DM halo. The 
mass distribution is then given by 



p(r, 9) = M.S{r) + Tv(r, 9) + p DM (r) 



(1) 



where M. is the black hole mass, T is the mass-to-light ratio, 
v is the stellar luminosity distribution, and pdm is the DM 
density profile. The potential <t>(r, 9) is derived by integrating 
Poisson's equation. 

In this potential, a representative orbit library is constructed 
that samples the phase space systematically. The genera- 
tion of the o rbit library i s descri bed in detail in lThomas et al.1 
(12004 and ISiopis et all (120091) . For the comparison with 
the data, we use a spatial grid of N r = 20 radial bins and 
Ne — 5 angular bins and use iV v = 15 velocity bins for the 
LOSVD at each spatial gridpoint. The galaxy potential and 
the forces are evaluated on a grid with 16 times finer reso- 
lution. For our axisymmetric code, there are three integrals 
of motion that sample the phase space accordingly: the en- 
ergy E, the angular momentum L : , and a non-classical third 
integral 73. The (E, L 7 )-plane is sampled based on the spa- 
tial binning dRichstone & Tremaindll98"8l) . We choose E and 
L z such that the respective orbits have their pericenter and 
apocenter in every pair of the ra dial grid bins. The thi rd inte- 
gral 73 is sampled as outlined by IThomas et al.l (120041) . First, 
orbits are dropped from the zero-velocity curve (defined by 
E = L 2 /(2r 2 cos 2 9) + O(r,0)), as in G03. This is done by 
using the intersections of the angular rays of the spatial grid 
with the zero-velocity curve as starting points for the integra- 
tion of the orbit's motion. However, this does not ensure a 
representative sampling of orbits. Such a sampling is indi- 
cated by a homogeneous coverage of the surface of section, 
i.e., the position of radii and radial velocities of orbits dur- 
ing their upward crossing of the equatorial plane. Therefore, 
additional orbits are launched to give such a homogeneous 
coverage. This method provides a complete sampling of the 
surface of section for given E and L- and thus a proper cover- 
age of phase space. We typically have 13,000-16,000 orbits 
in our library. The allocation of the individual orbits to the 
spatial grid points is based on the time they spend there. 

Given this orbit library, the orbit weights are chosen by 
matching the orbit superposition to the observed light distri- 
bution and the LOSVDs of the galaxy on the spatial grid. To 
fit the orbit li brary to the data, we us e the m aximum entropy 
technique of iRichstone & Tremainel ([1988). This method 
maximizes the function 

S = S - a X 2 , (2) 

where x 1 is the sum of the squared residuals over all spatial 
and velocity bins, e.g., 
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(3) 



Effect of a DM halo on the determination of dynamical black hole masses 



3 



Table 1 

Results for the Galaxy Sample 



Galaxy 


D (Mpc) 


M..G03 (Me) 


M/LG03 


A?.,noDM (M e ) 


Mj L uo om 


M.£im [M 


ol 


M/Ldm 


V c (kms- 1 ) 


r c (kpc) 


«inf (") 


(1) 


(2) 


(3) 




(4) 


(5) 




(6) 


V) 




(8) 


(9) 


(10) 


(11) 


JNOC oil 


25.5 


9.9 ±4.1 x 


10 7 


f~ Q 
0.0 


1.1 ±0.4 X 


10 s 


7.7 ± 0.5 


1.0 ± U.o X 


lu 


6.2 ± 0.7 


450 


14.0 


0.14 


JNOC Ilia 


24.2 


1.6 ± 1.0 x 


LO 7 


7.2 


1.4 ± 1.1 x 


10 7 


1 1.9 ± 1.1 


1.5 ± 1.5 X 


1 U 


11.8 ± 1.2 


300 


5.0 


0.02 


NGC 3377 


11.7 


1.1 ±0.6x 


10 s 


2.7 


1.6 ± 1.0 x 


10 8 


2.6 ±0.5 


1.9 ± 1.0 x 


10 s 


2.3 ±0.4 


350 


6.0 


0.69 


NGC 3384 


11.7 


1.8 ± 0.2 x 


10 7 


2.5 


8.0 ± 4.2 x 


10 6 


2.4 ±0.1 


1.1 ±0.5 x 


10 7 


2.2 ±0.1 


400 


8.0 


0.04 


NGC 3608 


23.0 


1.9 ± 0.9 x 


L0 8 


3.7 


4.6 ± 0.9 x 


10 s 


3.5 ±0.3 


4.7 ± 1.0 x 


10 s 


3.3 ± 0.3 


400 


10.0 


0.55 


NGC 4291 


25.0 


3.2 ± 1.6 x 


10 s 


6.0 


9.7 ± 2.0 x 


10 8 


6.0 ±0.5 


9.2 ± 2.9 x 


10 s 


6.0 ±0.8 


400 


8.5 


0.56 


NGC 4473 


17.0 


1.3 ±0.7 x 


10 8 


5.1 


5.9 ± 5.0 x 


10 7 


7.4 ±0.2 


1.0 ± 0.5 x 


10 s 


6.8 ±0.3 


400 


10.0 


0.15 


NGC 4564 


17.0 


6.9 ± 0.7 x 


10 7 


1.6 


9.8 ± 2.3 x 


10 7 


1.6 ±0.1 


9.4 ± 2.6 x 


10 7 


1.5 ±0.1 


350 


7.0 


0.19 


NGC 4649 


16.5 


2.1 ±0.6x 


10 9 


8.8 


3.9 ± 1.0 x 


10 9 


8.6 ±0.6 


4.2 ± 1.0 x 


10 9 


8.0 ±0.7 


500 


20.0 


1.51 


NGC 4697 


12.4 


2.0 ± 0.2 x 


10 s 


4.2 


2.2 ± 0.3 x 


10 s 


4.5 ±0.3 


2.0 ± 0.5 x 


10 s 


4.3 ±0.3 


450 


12.0 


0.45 


NGC 5845 


28.7 


2.9 ± 1.1 x 


10 s 


4.5 


4.5 ± 1.2 x 


10 s 


5.4 ±0.2 


5.4 ± 1.7 x 


10 s 


5.1 ±0.2 


300 


5.0 


0.30 


NGC 7457 


14.0 


4.1 ± 1.4 x 


LO 6 


2.8 


7.4 ± 4.2 x 


10 6 


2.7 ±0.5 


1.0 ± 0.6 x 


10 7 


2.6 ±0.5 


300 


5.0 


0.14 



Note. 



black hole mass and mass-to-light ratio using updated code but without including a DM halo in the dynamical modeling. Columns 7-8: black hole mass and 
mass-to-light ratio when a DM halo is included. Column 9-10: parameters of the circular velocity and the core radius for the logarithmic DM density profile 
used. Column 1 1 : radius of the black hole's sphere of influence, based on the black hole mass including the contribution of DM. 



where is the light in the kth bin, with the bins composed 
of the spatial position on the sky and the line-of-sight veloc- 
ity, thus the bin in the LOSVD, at that position. Hence, Ic is 
varying from 1 to Nd = N r NgN v . While /dat.A- refers to the 
measured light at that position, / mo d,t is given by the weighted 
sum of the contribution of all orbits to the fcth bin. 
S is the Boltzmann entropy 



5 =I> lo s(i)' 



(4) 



where w, is the weight of the individual orbit and V,- is the 
phase space volume of this orbit, i.e., the volume of the region 
in phase space that is represented by this orbit i, given by 



J": 



V = AEAL- T(r,v r )drdv, 



(5) 



where T(r, v r ) is the time between two successive crossings of 
the equatorial plane and AE and AL Z are the rang es in energies 
and angular momenta of the respective orbits dBinney et al.l 
fl98l lTTiomareTaDl200l . 

The parameter a in Equation (f2]i controls the relative weight 
of entropy and^ 2 for the maximization. We start with a small 
a, being dominated by the entropy maximization in the fit, 
and then iteratively increase it until there is no longer an im- 
provement in the x 2 - 

This procedure provides a value for the x 2 f° r on e com- 
bination of M., Y and DM halo. The best fit is found by 
the global minimum of x 1 for the variation of these param- 
eters. For the estimation of the parameter uncerta inties we 
adopt the usual Ax 1 = 1 criterion dPress et aD 19921) to obtain 
the 68% confidence intervals for one degree of freedom, thus 
when marginalizing over the other free parameters. 

4. RESULTS 

4.1. Models without a DM halo 

We first ran a set of models without including the contri- 
bution of a DM halo; thus, we set pdm = in Equation (|T). 
This assumption is consistent with most previous studies on 
black hole masses using dynamical models as well as with 
G03. As we are using the same data asthose of G03, the main 
difference is the improved modeling code. Thus, we would 



expect to recover similar black hole masses as in G03. We 
als o use slightly different distances to the galaxies, as given 
bv lGultekin etall (l2009bl) . 

For each galaxy we ran models on a fine grid in M. and 
mass-to-light ratio (M/L). Each model gives a best-fitting 
orbit superposition, and thus orbital structure for the galaxy, 
with a corresponding value of x 2 - The best fit M. and M/L 
are determined by the global minimum of the x 2 distribution. 

The x 2 distribution as a function of M. (marginalized over 
M/L) is shown as the blue line in FigureQ] We determined our 
stated best fit M. and M/L from their marginalized^ 2 distri- 
butions, using the mid-point of the Ax 2 = 1 interval, which 
corresponds to a lcr uncertainty. The results are presented in 
Table □ 



4.2. Models with a DM halo 

It has been shown that dynamical models are clearly 
able to detect and constrain the presence of a DM halo , 
if the data range suffic i ently far in radius j Rix et al] 119971: 
Kronawitter et al.ll2000t IThomas et al"1l2007t IWeijmans et alj 



20091; iForestell & Ge bhardt 2010). Due to the faintness of 



the stellar component at large radii, other kinematic tracers 
such as globular clu sters (Zepf et al. 2000; P ierce et al.l|2006 



Bridg es et al.l2006l) or pla n etary neb ulae (Mende z et al.12001 
Romanowsky et all 12003b ICoccato et alj 120091) have to be 



used. Furthermore, if dynamical coverage of both the central 
and the outer regions of the galaxy is present, it is possible to 
constrain M„ M/L, and the dark halo parameters by dynam- 
ical m odeling dGebhardt & Thomas! 120091 IShen & Gebhardt) 
I20l0h . 

Two common parameterizations for the DM halo are a 
Navarro-Frenk- White (NFW) profile (iNavarro et al.l 119961) 
and a DM distribution based on a cored logarithmi c poten- 
tial dBinnev & Tremainel [19871 IThomas etall 12005b. For 

j.lTho " 



sample of 17 early-type galaxies. IThomas et al. (2007) found 
both profiles to give consistent re sults, with tentative evidence 
to favor a logarithmic dark halo. iGebhardt & Thomas! (120091) 
confirmed this result for M87, and iMcConnell et al.l ( 120101) 
found consistent results for M. using either a logarithmic dark 
halo or an NFW profile. In the following, we will use a DM 
halo with a cored logarithmic potential, whose density profile 
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M. [M e ; xlO 8 M, [M Q ] xlO 7 M. [M Q ] xlO 1 




0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.5 1.0 1.5 2.0 

M, [M Q ] xlO 7 M, [M Q ] xlO 9 M. [M Q ] xlO 9 





1 2 3 4 5 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 2.5 

M. [M Q ] xlO 8 M. [M ] xlO 9 M, [M ] xlO 7 



Figure 1. Comparison of the x 2 distributions as a function of M, (marginalized over MIL). The black solid line shows the models including a DM halo. Their 
X 1 values always show the actual modeling result. The \ 2 distribution for the models without a DM halo is shown as the blue dashed-dotted line. The zero point 
has been shifted for NGC 2778 and NGC 3384 by an offset given in the figure as A^ 2 (in blue). The^ 2 distribution of G03 is shown as the red dashed line, offset 
by the value given as AXq 03 (in red). The distributions have been scaled in M., to account for the difference in the assumed distance. 
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Figure 2. Two-dimensional contour plot of x 2 as a function of M. and Ml L for each galaxy. Here the models with a DM halo included are shown. The contour 
lines show Ax' values of 1.0, 2.71, 4.0, and 6.63 (corresponding to 68%, 90%, 95%, and 99% confidence for one degree of freedom). The points indicate the 
individual models we ran. The contours are derived from a smoothing spline to these models. The cross represents the best fit. 
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is given by 



PdmO) 



V 



3r 2 + r 2 



4ttG (r 2 + r 2 ) 2 



(6) 



where V c is the asymptotically constant circular velocity and 
r c is the core radius, within which the DM density is approxi- 
mately constant. 

Our data in general do not constrain the DM profile, as we 
are lacking kinematic information at large radii. While for a 
few galaxies in our sample, large radii kinematic information 
for the stars, globular clusters or planetary nebulae exist in 
the literature, we do not include them in this analysis. In this 
work, we are not aiming at constraining the DM halo itself, 
but we are mainly interested in the effect of including such a 
halo for the recovered black hole mass. We leave a more de- 
tailed investigation of the combined DM halo and black hole 
properties for these individual galaxies to future work. This 
also allows a better direct comparison to the work of G03 and 
the models without a DM halo, presented in the previous sec- 
tion. 

Therefore, we assume a fixed DM halo, with fixed param- 
eters V c and r n . These are taken from the scaling relations 
presented bv lThomas et al.l d2009t) . based on the galaxy lumi- 
nosity: 

log r c = 1.54 + O.63(log(L B /L ) - 1 1) (7) 

log V c = 2.78 + 0.21(log(L fl /Zo) - 11) , (8) 

and given in Table Q] These scaling relations have been es- 
tablished based on a sample of 12 early-type galaxies in the 
Coma cluster with old stellar populations. While our sam- 
ple does not have to follow these scaling relationships ex- 
actly, they at least provide well motivated parameters for the 
DM halo. Younger early-type galaxies and disk galaxies have 
been found to have on average a less massive halo, thus our 
approach tends to maximize the DM contribution. We inves- 
tigate the effect of changing the assumed DM halo on the cen- 
tral black hole mass further below. 

Thus, for each galaxy, we ran a grid of models for vary- 
ing M. and M/L with fixed DM halo. We show the two- 
dimensional distribution of x 2 as a function of M. and Mj L 
in Figure [2] The contours are based on the x 2 values of the 
underlying grid points, ap plying a two-dimensional smooth- 
ing spline (lDierckxlll993l) . The marginalized^ 2 distribution 
as a function of M. is shown as the black line in Figure Q] We 
again determined the best-fit values for M, and M/L from the 
marginalized distribution and have given them in Table Q] 

5. COMPARISON OF BLACK HOLE MASSES 

5.1. Comparison with \Gebhardt et al.\ f feOOil) 

As we are using the same data as those in G0 3, th e only 
difference between the work presented in Section 14.11 and in 
G03 is the improved modeling code. Thus, we would expect 
to recover the same black hole masses as in G03. 

For a comparison with G03, their masses are fi rst increased 
by a f actor of 1 .09, due to a unit conversion error CSio pis et al.l 
2009), and then are rescaled, according to the difference in 
the ado pted distance, a s suming M. cc d. These masses are 
listed in iGiiltekin et al.l (l2009bl) . apart from NGC 821. This 
galaxy has an erroneou s black hole mass in G03, corrected in 
iRichstone et alj ( 120041) . After accounting for the factor of 1 .09 
and the distance difference, the black hole mass for NGC 821 
is M. = 9.9 x 1O 7 M . 
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Figure 3. Black hole masses given by Gebhardt et al. 12003 ) vs. the black 
hole masses determined in this work (without a dark halo). The solid line 
shows a one-to-one correspondence. 
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Fi gure 4. Ratio between the M. determined in this work and the M. given 
m rGebhardt et aT] <2003l) , as a function of the G03 M.. The blue squares 
are for the models without a dark halo, while the black circles correspond to 
the result, when a dark halo is included in the models. The solid line is a 
one-to-one correspondence to the G03 values. 

In Figure[3]we compare the black hole masses, determined 
without including a DM halo, with the black hole masses 
given in G03. In Figure [4] we show as blue squares the ra- 
tio between both mass determinations as a function of the G03 
mass. The marginalized ;f 2 distributions for the individual ob- 
jects are shown in FigureQ] as blue dashed-dotted lines for the 
current work masses and as red dashed lines for the G03 dis- 
tributions. Note that the G03 distributions are offset in x 2 , 
such that the minimum corresponds to the minimum of the^ 2 
distribution including a dark halo, shown in black. The reason 
for the offset in^ 2 is mainly due to the larger number of orbits 
used in the current modeling, compared to G03. 
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For three objects (NGC 821, NGC 2778, and NGC 4697) 
the difference in M, is less than 20 %, thus consistent with our 
previous work. The internal structures of the dynamical mod- 
els (as discussed below) are similar for these three galaxies in 
the old and new models, which explains the reason for lack of 
change. The small difference is probably due to the presence 
of numerical noise in the models. This noise is mainly caused 
by the use of a finite number of individual orbits instead of 
a smooth orbit distribution function. The comparison of the 
X 1 distribution for the three objects shows that they are basi- 
cally consistent, while the distribution may widen, possibly 
due to a more complete orbit library. Also, for NGC 3377 and 
NGC 4564 ,the difference in M, is within the stated uncertain- 
ties. 

However, for the rest of the objects the new M. is sig- 
nificantly offset from the previous estimate, not simply ex- 
plained by numerical noise. For two objects, NGC 3384 and 
NGC 4473, M. decreases; for the rest there is an increase in 
M„ by up to a factor of three. The mean increase for this 
sample is a factor of 1.46 with the standard deviation of 0.73. 
This r esult is in line with the findings of IShen & Gebhardtl 
(120101) . who found an increase in the mass by a factor of two 
for NGC 4649. This o bject is included in Figur es [3] and Has 
the highest mass point. IShen & Gebhardtl (120101) attribute the 
mass difference to the better orbit sampling in the new mod- 
els. In particular, they argue that the old orbit sampling lacks 
high energy, nearly circular orbits, which lead to an underes- 
timate of M,. NGC 4649 is a core galaxy, and it is important 
to note that all galaxies with an increase in mass by more than 
a factor of two are core galaxies as well. This seems to indi- 
cate that the previous orbit sampling was not able to properly 
model core galaxies. 

To investigate this issue further for our whole sample, we 
inspect the internal orbit structure, looking for any clear dif- 
ference between the models. To do so, we examine the shape 
of the velocity dispersion tensor, represented by the ratio of 
radial to tangential dispersion cr r lcr t . The tangential disper- 
sion includes contributions from random as well as from or- 
dered motion; thus, it is given by crj = cr 2 g + cr 2 + V?. In 
Figure [5] we compare the internal dispersion ratio cr r lcr t for 
the best-fit models presented here, with and without a DM 
halo (as blue solid and black dashed dotted lines, respec- 
tively), with the ratio for the models in G03, shown as red 
dashed lines. We also indicate the black hole sphere of in- 
fluence R m f = GM.cr~ 2 , assuming the new M, (without DM 
halo). The galaxies with consistent black hole masses, such 
as NGC 3377 and NGC 4697, also exhibit consistent inter- 
nal structure. On the other hand, the galaxies with the largest 
mass increase, especially the core galaxies such as NGC 4291 
and NGC 3608 show a clear difference in the internal struc- 
ture. First, there is a strong radial bias at large radii for these 
galaxies, especially compared to the previous dispersion ra- 
tio. However, this radial bias is mainly outside the range for 
which kinematic data are available and is therefore driven by 
the maximization of the entropy. We do not expect these or- 
bits to have an influence on the black hole mass determination. 
Second, there is a stronger tangential bias inside the black 
hole sphere of influence. In particular, the previous models 
exhibit a radial bias within R m f for the largest outliers. G03 
only sampled the zero-velocity curve, instead of the whole 
phase space, and due to the coarse sampling of drop points 
they missed the orbits near the pole that are nearly circular. 
This sampling then causes a radial orbital bias. This radial 



bias is removed in the models presented here, using a better 
orbit sampling. An increase in tangential orbits will reduce 
the projected line-of-sight velocity dispersion and therefore a 
more massive black hole is required to match the observed ve- 
locity dispersion profile. On the other hand, NGC 4473 which 
shows a decrease in the determined black hole mass, exhibits 
a stronger radial bias in the new modeling compared to G03. 
This radial bias is probably caused by the presence of a nu- 
clear disk in this galaxy. 

Thus, we find that the main reason for the change in black 
hole mas s is the different orbit s ampling used, as already 
found by IShen & Gebhardtl J2010h for NGC 4649. We now 
cover the phase space more completely and therefore also in- 
clude orbits missed by the previous sampling. This issue is of 
special importance for core galaxies, as they often show a sig- 
nificant tangential orbital bias in their center, i.e., they usually 
have the largest <x, (G03). 

To illustrate this point, we computed the difference of the 
dispersion ratio between G03 and this new model in a shell 
inside the black hole sphere of influence: 



ARo 



£ 



[(07/(77)003 - (ov/cr,)] dr . 



(9) 



with fmin = 0.1/?i n f and r max = R m f. This quantity is just a 
simple and quick way of quantifying the change in the orbital 
structure and is just meant to highlight the relation between 
the change in orbital structure and the change in black hole 
mass. In Figure [6] we plot it against the ratio of the black 
hole masses M./M.G03. There is a clear correlation between 
the quantities, confirming our previous argument. We have 
also tested the effect of decreasing the number of orbits in the 
modeling, but saw no clear influence on the best-fit black hole 
mass. Thus we confirm our previous results in finding that the 
reco vered black hole mass is not affected by the number of or- 
bits (|Gebhardtll2004lRlchstone et all2004t IShen & Gebhardtl 
I2010l 



This investigation emphasizes the need for a complete or- 
bital sampling of phase space for dynamical modeling of 
galaxies, especially of core galaxies. 

5.2. Effect of a dark matter halo on the determined black 
hole mass 

The main motivation of this paper is to investigate the ef- 
fect of the inclusion of a DM halo in the dynamical modeling 
on the determined black hole mass. In Figure|7]we show the 
difference in black hole mass between the models with and 
without the inclusion of a DM halo as a function of the res- 
olution of the black hole sphere of influence divided by the 
spatial resolution of the kinematic observation {R m ild KS ). For 
the computation of R m f we used the black hole mass including 
a DM halo (given in Table [TJ. The spatial resolution is given 
by the seeing and the aperture of the HST kinematic observa- 
tions; thus, d res = 0.08 for the STIS data and d KS = 0.15 for 
the FOS data. 

As expected, there is a general trend of an increase in M, 
when a DM halo is included. For five objects, we find almost 
no change in M., while for one object - NGC 2778 - the sig- 
nificance of the black hole detection even vanishes, with the 
minimum^ 2 for no black hole. The other six galaxies show 
an increase in the measured M, between 20% and 80% when 
a DM halo is included. The most extreme case is NGC 4473 
probably due to the presence of a nuclear disk, with an in- 
crease of a factor of 1.8 when a DM halo is included. For the 
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Figure 5. Ratio of the internal velocity dispersions as a fu nction of radius for the best model of each galaxy. Shown is the ratio <T r jcr t along the major axis. The 
result with a DM halo, without a DM halo, and the result of|Gebhardt et al. 1 2003) (without a dark halo) are shown as the black dashed, blue solid and red dashed 
lines, respectively. The horizontal dashed line corresponds to a non-rotating isotropic model. The vertical dashed-dotted line indicates the black hole's sphere of 
influence, assuming the black hole mass, determined without including a DM halo. The vertical dashed lines show the radial extent of the ground-based data. 
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Figure 6. Ratio betwee n the M. determined in this work and the M. given 
in lGebha rdt et al. 1 2003), as a fu nction of the excess in radial motion in the 
internal structure of the models of Gebhardt et al. 1 2003), as defined by Equa- 
tion {5J. A correlation between both quantities is apparent. 

whole sample we find a mean increase of a factor of 1 .22 with 
standard variation of 0.27. This increase is much less than the 
factor of more than two found for M87. In contrast to M87, 
our data set contains no stellar kinematic information at large 
radii but includes HST data at small radii. Thus, we are better 
able to probe the region affected by the presence of the black 
hole at the center. 

In Figure [7] there appears to be a trend of a larger bias for 
objects where R m f is less well resolved, as would be expected. 
However, due to the black hole mass uncertainties there is 
no statistically significant relation. The most massive galaxy 
in our sample, NGC 4649, is not shown in the figure, as it 
would appear at Ri n f/d les a 20 with no significant change 
in black hole mass. M87 would lie at Ri n f/d KS » 1.5 and 
W..dm/-M.,iioDM ~ 2.8. In contrast to M87, the galaxies in our 
sample with a less well resolved sphere of influence exhibit- 
ing a smaller change in the determined M. are less massive 
and probably reside in less massive DM halos. This indicates 
that especially for massive galaxies properly resolving /?j n f is 
important to determine M. under the consideration of DM. 

In Figure Q] the marginalized x 1 distributions for the in- 
dividual objects with (solid black line) and without (dotted 
dashed blue line) a DM halo are shown. For the five ob- 
jects with almost no change in M. (NGC 3377, NGC 3608, 
NGC 4291, NGC 4649, and NGC 4697), there is also no 
change in the y 1 (apart from NGC 4649). For the other 
galaxies, including a reasonable DM halo improves the fit in 
terms of^ 2 . The most convincing cases are NGC 821 and 
NGC 3384, where the model without a DM halo is excluded 
at more than 3<x significance. Thus, while we are not able to 
constrain the shape of the DM halo, at least for some galax- 
ies the presence of such a halo is supported. In total, for 
six galaxies (NGC 821, NGC 2778, NGC 3384, NGC 4473, 
NGC 4564, and NGC 5845), the model without a DM halo is 
excluded with at least 2<x significance. 

The mean increase of the black hole mass goes along with a 
decrease of the mass-to-light ratio, as expected. This indicates 
the degeneracy present between the stellar mass-to-light ratio 
and the DM contribution in dynamical models. For the whole 
sample we find a decrease in M/L of 6% with a scatter of 5%. 

Even if our choice of DM halo is well motivated by the 




Figure 7. Ratio between the M. without accounting for the DM halo and the 
M. including a DM halo in the model, as a function of the ratio of black hole 
sphere of influence over the spatial resolution. The arrow indicates the upper 
limit for the mass of NGC 2778 when DM is included. The solid line shows 
a one-to-one correspondence between both masses. 



scaling relations of iThomas et al.l (|2009), it is basically an ad 
hoc assumption we had to make as we do not have the data 
to robustly constrain the DM halo profile. To at least test the 
effect of changing the assumed DM halo on the black hole 
mass, we ran a set of models, changing V c in the logarith- 
mic DM potential. We restrict ourselves to changing only this 
one parameter, as we want to avoid sampling the whole four- 
dimensional parameter space. It has also been found that V c 
and r c a re degenerate, especially if the large rad ii coverage 
is poor dShen & GebhardtH20ld: [Forestelll 120091) . For each 
galaxy, we assume a twice as massive DM halo and a DM 
halo about half as massive, as well as some additional values. 
The results are shown in Figures [8] and [9] We confirm the 
basic trends of an improved^ 2 for a reasonable massive halo 
and an increase in M, for a more massive halo. The range 
of given M, approximately covers the range consistent with 
the current data, as long as the DM halo is not constrained for 
these galaxies. 

5.3. Notes on individual galaxies 

In the following, we provide more detailed information on 
the black holes for some individual galaxies. 

NGC 821. There are two recent stu dies on the DM halo of 
this galaxy, providing large radii data. IWeijmans et al.l (120091) 
used SAURON data to measure LOSVDs out to ~ 4 effec- 
tive radii. Their assumed DM halo gives M DM = 9x 1O 9 M 
within the eff ective radius (assuming R,, = 5.1 kpc), using an 
NFW profile. iForestell & Gebhardtl OOlOh used long-slit data 
from the Hobby-Eberly Telescope to measure the LOSVD out 
to ~ 2R e . Assuming their power law fit to the DM halo, we 
find M DM = 8 x 10 9 M X within R e . Our assumed DM halo 
is more than twice as massive. Thus, M, for the true DM 
halo should be contained within the range spanned by our no 
DM and DM solution. Including these large radii data into 
the dynamical models is beyond the scope of this paper. The 
nuclear supermassive black hole in NGC 821 has been de- 
tected as a weak X-ray sour ce, implying a very w eak level 
of activity (L x /L Edd ~ 10' 8 ; lPellegrini et al.ll2007h. There is 
also evidence for the presence of a iet dFabbiano et alj|2004l 
iPellegrini etaT]|2007l) . 
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Figure 8. Variation of the black hole mass and the minimal x 2 as a function of the assumed DM halo for each galaxy. In the left columns of the subpanels the 
change in M. for various values of Vdm is shown. The error bars correspond to lcr (i.e., A^ 2 = 1). In the right columns the corresponding value of x 1 is shown. 
A spline interpolation is shown as the solid line. 



NGC 2778. This galaxy already had the least confident 
black hole detection in G03. Assuming the value for M. of 
our no DM halo model, we do not resolve the black hole's 
sphere of influence (Ri n f/d les = 0.2). Including a DM halo in 
the model improves the fit significantly, but the significance 
of the black hole detection disappears. However, the previous 
M. estimate is still fully consistent with the lcr upper limit of 
M., up = 2.99 x 1O 7 M that we derive for NGC 2778 under 
the presence of a DM halo. This behavior might indicate the 
need to properly resolve 7?; n f when a DM halo is included to 
properly determine M.. 

NGC 3377. The black hole mass for this galaxy increased 
by ~ 70%, compared to G03, mainly caused by the stronger 
widening of the confidence contours at the high-mass end than 
at the low-mas s end. The previou s value is still fully consis- 
tent within 1 cr. ICopin etaD (120041) reported a black hole mass 
of M. — 8.3 x 10 7 M (for our assumed distance) based on In- 
tegral Field Unit (IFU) observations with SAURON and OA- 
SIS, also still consistent with our results within lcr. The first 
detection of a black ho le in NGC 3377 has been reported by 
iKormendv et al.l d 19981) . based on ground-based observations. 
Using an isotropic model, they found M. = 2.1 x 10 8 M and 
M/Ly = 2.0 (for our assumed distance), in good agreement 
with our results. NGC 3377 is a rapid rotator and close to 
isotropy, justifying the isotropic assumption for this galaxy. 
NGC 3 377 exhibits a nuc lear X-ray source, showing a jet like 
feature (ISoria et al.ll2006l) . 



NGC 3384. Besides NGC 2778, this is the only other galaxy 
for which the sphere of influence is not resolved. While for 
NGC 2778 the new code does not lead to a change of the x 1 
distribution, for NGC 3384 M. decreases. This might indicate 
a larger uncertainty in the determination of M, using different 
modeling codes when R m { is not resolved. NGC 3384 is the 
galaxy with the strongest constraints on the presence of a DM 
halo. For this galaxy, we ran a grid of models changing r c as 
well as V c , but we found no change in^ 2 for different values 
of r c . However, we are able to set a lower limit on V c with 
V c >~ 350 km s -1 at lcr confidence. The no DM halo model 
is excluded at more than 3cr confidence (see Figure©. 

NGC 4473. NGC 4473 shows evidence for a central stellar 
disk both in the imaging and the kinematics, as discussed by 
G03. We followed G03 and include a central exponential disk 
and also assumed a galaxy inclination of 71°, as found for the 
disk component. Thus this galaxy is the only case in our sam- 
ple not modeled with an edge-on inclination. The presence of 
the disk has a distinct influence on the measured black hole 
mass causing a relatively large difference between the models 
with and without a DM halo. 

NGC 4564. This galax y is known to have a nuclear X-ray 
source (ISoria et al J 120 06), indicating the presence of an ex- 
tremely sub-Eddington accreting AGN. 

NGC 4649. This ob ject has recently been studied by 
IShen & Gebhardt! (120101) including a DM halo in the models. 
In addition to the stellar kinematics used in this work, they 
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Figure 9. Same as Figure [8]for the second half of the galaxy sample. 

included globular cluster velocities from H wang et al.l (120081) . 
Thus, our results are not directly comparable. They report 
values of M. = (4.5 ± 1.0) x 10 9 when including a DM halo 
in the models and M. = (4.3 + 0.7) x 10 9 without a DM halo. 
Our results for M, are consistent with their work. 

NGC 4697. The black hole mass for this galaxy is basically 
unchanged using the modifie d code and includ ing a DM halo. 
This result is consistent with iForestelll ( 120091) . She used the 
same data and model code as we did, but augmented by kine- 
matics of pl anetary nebulae at large radii dMendez et al.l200"Tl 
120081 12009b . constraining M. and the DM halo at the same 
time. Her best-fit model has M. = 2.1 x 10 s M , M/L = 4.35, 
V c = 388 km s _1 , and r c — 9 kpc, assuming a logarithmic 
halo. We find identical values for M. and M/L, using slightly 
different DM halo paramet ers. NGC 4697 h as a nuclear point 
source detected in X-rays (^Soria et al. 2006), showing that its 
black hole is active at a low rate. 

NGC 5845. While there is a moderate increase in M, when 
a DM halo is included in the model, the^- 2 distribution flattens 
at the high-mass end, due to an increased degeneracy between 
M. and M/L. Increasing the mass of the DM halo strongly 
enhances this degeneracy, leading to an almost unconstrained 
M. over a wide mass range, until for V c ~ 600 km s _1 the min- 
imum switches to M, « 1 .7 x 10 9 M , still with strong degen- 
eracy between M, and M/L. Kinematic data at large radii, to 
better constrain the DM halo and M/L would be desirable for 
this galaxy. The model without a DM halo is excluded with 
more than 3cr significance for NGC 5845. There is a nuclear 
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X-ray source here as well dSoria et al.l 120061) . There is evi- 
dence for obscuration of the black hole by a dusty disk, with 
the X-ray emission originating from scattering of the AGN 
continuum emission on the surrounding plasma. 

6. THE BLACK HOLE-BULGE RELATIONS 

As our sample constitutes a significant fraction of the 
galaxy sample for which dynamical black hole masses are 
available, it is worth looking at the effect of these new black 
hole mass measurements on the black hole-spheroid relations, 
namely the M. - cr t and M. - Ly relationships. We used 
the sample of |Gultekin et al. ( 2009b) as the reference sample, 
containing 49 M. measurements and 18 upper limits, includ- 
ing our 12 objects. 

For the fitting, we used a generalized maximum likeli- 
hood rnethodas described by [GUtekin et al. (2009b, ; see 
also IWoo et afl (120101) 1. We' minimize the likelihood func- 
tion S = -2 In .£, with £, = Y\i liiMi, si) being the product of 
the likelihoods for the individual measurements of black hole 
mass yu = log M, and bulge property s = log cr« or s — log Ly. 
The likelihood for measuring the mass /i, and bulge property 
Si for given true mass fi and true bulge property s is: 

Si ) = J | n)Q s ( Si | s)PQi | s) dn ds . (10) 

We assume Q M , Q s and P to have a log-normal form, with ctq^ 
and ctq corresponding to the measurement uncertainty in the 
black hole mass and bulge property, and <x P = eo is the intrin- 
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sic scatter in the black hole mass-bulge propert y relation. Up- 
per limi t; 

(2009b). Thus, we minimize 



?er limi ts are incorporated in the fit, following Giilte kin et al.l 
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Figure 10. Left panel: M. - <x„ relationship. The red squares show our sample with the M. values given by[Gebhardt et al. 1 2003), the black circles give the M. 
determine d in this work with the i nclusion of a DM halo. The black solid line shows our updated result for the M. - cr, relation, while the dashed line shows the 
relation byJGtiltekin et al. (2009b). The crosses show the rest of their sample. Right panel: M. - Ly relationship. The symbols are the same a s in the left panel. 
The solid line is our best-fit when including a DM halo for our 1 2 galaxies. The dashed line and crosses are again from lGultekin et alj <2009bD . 

malization increased slightly as well as the i ntrinsic scatter 
in both relations. We also fitted the sample of iGiiltekin et alj 
d2009bl) with our updated black hole masses, without account- 
ing for a DM halo. Most of the change in the M, - cr, and 
M. - Ly relationships is caused by the improved masses. The 
effect of the inclusion of a DM halo on these relationships is 
marginal. 

However, this is not a full correction of the black hole- 
bulge relationships for the effect of a DM halo on the black 
hole masses, as it is restricted to our sample of 12 galax- 
ies. The rest of the galaxies with stellar dynamical black hole 
mass measurements potentially suffer from the same system- 
atic bias. Ideally, a correction would consist of a re-modeling 
of these galaxies including a DM halo, as performed in this 
work for the sample of G03. However, we can use Figure [7j 
as a guideline for an average correction. Figure [7j indicates 
that the correction factor depends on the resolution of the 
sphere of influence. For Ri n f/d KS > 3 including or ignoring 
a DM halo in the modeling gives consistent results, while for 
lower values there is on average a systematic bias with a mean 
( M, dm /M m _ no dm ) = 1-5 for our sample. 

To estimate the effect on the black hole-bulge relations, we 
increased all stell ar dynamical bl a ck hole mass measurements 
in the sample of IGiiltekin et all (l2009bl) with R m f/d KS < 3 
by this average factor and re-fitted the relations. To inves- 
tigate the pure change due to the DM halo, we also fitted 
the blac k hole-bulge relations to the sample of IGiiltekin et alj 
(l2009bl) . but with M. of the 12 galaxies of our work replaced 
by our results without a DM halo. Compared to the best-fit 
to this sample, we found a slightly increased slope, a con- 
sistent intrinsic scatter, and an increase in normalization by 
0.04 dex. We found a normalization, slope, and intrinsic scat- 
ter of (8.21,4.38,0.42) for the M. - rx* and (9.05,1.07,0.41) 
for the M, - Ly relationship. 

Additionally, we fitted only our sample with the values for 
M. with and without including a DM halo in the models. We 
recovered an increase in the normalization of ~ 0.07 dex, cor- 
responding to the mean increase in M, in the sample, while 
the slope is consistent and the intrinsic scatter decreases. 



with a and p being the normalization and the slope of the 
black hole-bulge relations, ef ot ; = ct 2 q i + ct 2 q i + e^, N is the 
number of black hole measurements, M is the number of up- 
per limits, and l n \.j is the likelihood of the upper limit as in 
IGiiltekin et all d2009bl 

We first fit the M. - cr t and M. - Ly relationships using 
the sample of IGiiltekin et al.ld2009bl) . finding identical results. 
We then updated their black hole masses with our new values 
for the 12 objects in our sample. We find 

log(M./M ) = (8.18+0.06) + (4.32+0.31) log(o-„/200 km s" 1 ) 

(12) 

with intrinsic scatter eo = 0.44 + 0.06 and 

log(M. /M ) = (9 .0 1 +0. 1 0) + ( 1 .06 +0. 1 5) log(L v / 1 1 1 L e v ) 

(13) 

with intrinsic scatter eo = 0.41 + 0.04. They are shown in 
Figure [TOl 

Note that bo th relations are not bas ed on exactly the same 
samples. As in Giiltekin et al. 1 2009b). for the determination 
of the M. - Ly relationship we restricted the sample to ellip- 
tical and SO galaxies with reliable bulge-disk decomposition. 
When using the same restricted subsample for the M. - cr t 
relationship, we find a shallower slope (J3 = 3.80 ± 0.33) and 
a reduced intrinsic scatter (eo = 0.34 + 0.05), lower than for 
the M. - Ly relationship for the same sample. Restricting the 
sample in this manner is supported by observations that sug- 
gest that s piral galaxies do no t follow the M, - cr t relation of 
ellipticals dGreene et al.ll2010l) . 

We also used a generalized least squares method to incor- 
porate measurement unce rtainties in both variab les and intrin- 
sic scatter as described in lTremaine et alj d2002l) . omitting the 
upper limits, which y ields consistent results. Compared to 
IGiiltekin et al.l d2009bl) . we find only a slight change for the 
best-fit. While the slope of the relation is consistent, the nor- 
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7. CONCLUSIONS 

We investigate the influence of accounting for the presence 
of a DM halo in the stellar dynamical modeling of galaxies 
on the measured black hole m asses. We use a samp le of 12 
galaxies, already analyzed by iGebhardt et al.l (120031) . which 
have ground based as well as high-resolution HST observa- 
tions of the stellar kinematics to address this issue. 

We model these galaxies without and with the presence of 
a DM halo. In the first case we found a significant differ- 
ence of the measured black hole masses between our previous 
results for a large fraction of the sample. For most of the 
objects the ma s s incr eased compared to the values given by 
IGebhardt et al.l (120031) . We ascribe this difference to the im- 
proved code, exhibiting a better coverage of the phase space 
for the generated orbit library. This shows the importance of 
a dense coverage of phase space in the dynamical models. 

Second, we include a reasonable DM halo into the models, 
using a scaling relat ionship based on the galaxy luminosity 
dThomas et al.l 120091) . We find an increase of the measured 
black hole mass, but much less than what has been found for 
M87 and NGC 6086. For these two galaxies, kinematic infor- 
mation is available only at large radii, whereas for our sample 
we have high-resolution data covering the central parts of the 
galaxies. Thus, the black hole mass is better constrained by 
central kinematic observations and less affected by the pres- 
ence of a DM halo in the models. 

Using different massive DM halos for the same galaxy, we 
confirm the trend of an increase of the recovered black hole 
mass for a more massive halo as well as a decrease of the 
mass-to-light ratio. Based on ax 2 analysis, the presence of a 
DM halo is implied for five of the 12 galaxies with at least 2cr 
significance, although we are not able to constrain the shape 
of the DM halo. 

We study the consequence of our new black hole mass mea- 
surements on the M. — q-„ and M. - Ly relationships, updat- 
ing the sample of Gultek irTet alj d2009b ) with our results. We 
found only a mild change in the best-fit values, still consis- 
tent with the previous estimate, with a slight increase in the 
normalization and the intrinsic scatter. We estimated the total 
effect of a black hole mass increase for galaxies studied by 
stellar dynamics by accounting for a DM halo which will lead 
to an increase in the normalization by ~ 0.04 - 0.07 dex. 

Even if our sample shows only a mild influence of the DM 
halo on the black hole mass, a DM halo is clearly present. 
Thus it is necessary to take it into account in the modeling of 
the galaxy to avoid a systematic bias. 
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